Entropy Driven Dimerization in a One-Dimensional Spin-Orbital Model 
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We study a new version of the one-dimensional spin-orbital model with spins S — 1 relevant 
to cubic vanadates. At small Hund's coupling Jh we discover dimerization in a pure electronic 
system solely due to a dynamical spin-orbital coupling. Above a critical value Jh, a uniform 
ferromagnetic state is stabilized at zero temperature. More surprisingly, we observe a temperature 
driven dimerization of the ferrochain, which occurs due to a large entropy released by dimer states. 
This dynamical dimerization seems to be the mechanism driving the peculiar intermediate phase of 
YV0 3 . 

PACS numbers: 75.10.Jm, 75.40.Mg, 05.70.-a 



In transition metal oxides with perovskite structure the 
d-orbital level splits into a triplet of ti g symmetry and an 
e g doublet at higher energy. The sign and magnitude of 
the spin-spin interactions is determined by the orbital 
occupation resulting in a strong spin-orbital coupling. 
This leads to unusual magnetic properties of many vana- 
dium, titanium and manganese oxides JjJ. The simplest 
spin-orbital model for spin-1/2 objects with twofold de- 
generated orbitals is the SU(4) symmetric model, where 
Hund's coupling, orbital anisotropy as well as the Jahn- 
Teller effect are ignored . The one-dimensional version 
of this model has been solved by Bethe ansatz || . The 
ground state and the excitation spectrum consisting of 
three gapless mixed spin-orbital modes is obtained. The 
idea of mixed spin-orbital excitations was also crucial 
to understand anomalous magnetic properties of LaTiOs 
a system showing a coherent orbital- liquid ground 
state H). Thermodynamic properties of the SU(4) sym- 
metric model in ID have been investigated by the quan- 
tum Monte Carlo method |J. 

Recently vanadium oxides, where the ti g orbitals 
are filled by two electrons, have attracted much atten- 
tion, because of their unusual magnetic properties in- 
cluding temperature-induced magnetization reversals in 
YVO3 @. YVO3 shows a very puzzling intermediate 
phase at 77 < T < 118K with C-type magnetic order 
(ferromagnetic chains along the c axis with antiferromag- 
netic coupling between the chains), which changes to a 
conventional G-type Neel structure at T = 77K through 
a first order transition || . A significance of ti g orbital de- 
grees of freedom for understanding the peculiar magnetic 
behavior of YVO3 has been emphasized [§ [lCj H) . 

The relevant spin-orbital model to describe magnetism 
of vanadium oxides has to involve S = 1 spins because 
Hund's coupling Jh is large [[l3| . For the cubic vanadates 
LaVC>3 and YVO3 a realistic superexchange model has 
been derived in Ref. |l2| where two out of the three ti g 
orbitals are active along a given cubic axis leading to a 
t = 1/2 orbital pseudospin. In the classical limit S ^> 
1 and neglecting Hund's coupling Shen et al. (Ml have 



shown that the ground state in ID is an orbital valence 
bond (OVB) solid formed by orbital singlets and parallel 
spins, where neighboring OVB's are noninteracting. 

In this paper we treat a realistic quantum spin-orbital 
model with S — 1 and a twofold orbital degeneracy in 
ID where the effects due to Hund's coupling are included. 
Thermodynamics of two distinct phases, a 4-site periodic 
quantum dimer phase, and a ferromagnetic state that oc- 
cur at larger Jh, is studied numerically. A critical value 
of Jh separating the two phases is obtained. The most 
important observation is that strong dimer correlations 
develop at finite temperature on the ferromagnetic side of 
the transition, which are driven by large entropy of low- 
lying dimer states of the model. We discuss the relevance 
of these findings to the intermediate phase of YVO3. 

The Hamiltonian of the model is given by 



(1) 



where S being an S — 1 spin operator, and J = 4t 2 /U 
represents an overall superexchange energy scale. The 
operators Jij and Kij describe orbital exchange pro- 
cesses on the bond 
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with r acting in a r = 1/2 orbital pseudospin space. 
The coefficients R = 77/(1 — 3?y) and r = 77/(1 + 277) 
originate from the Hund's rule splitting of the excited 
rj 2g multiplet by 77 = Jh/U. This is the Hamiltonian 
derived for cubic vanadates in Ref. ]l2| restricted here 
to one cubic axis. Physically, pseudospin r describes ID 
fluctuations of xz and yz orbitals along the c axis, while 
the third orbital xy is frozen below a structural transition 
at ~ 200K g, 0. 

To study the thermodynamic properties we apply the 
transfer-matrix DMRG, a numerical method which is 
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particularly suited because the thermodynamic limit is 
performed exactly. By means of the Trotter formula the 
partition function is expressed as a product of local trans- 
fer matrices evolving along the spatial direction. For the 
calculations we use the well known decomposition of the 
Hamiltonian into even and odd parts leading to a clas- 
sical model with checkerboard structure jl5| as well as 
a different Trotter-Suzuki mapping leading to a classical 
model with alternating rows |l6|]. Within the alterna- 
tive mapping correlation lengths (CLs) and correspond- 
ing wavevectors are determined unambiguously by the 
eigenvalues of the quantum transfer matrix (QTM). In 
the thermodynamic limit the free energy in both cases is 
given solely by the largest eigenvalue of the QTM, which 
is real and unique. The infinite DMRG algorithm is used 
to decrease the temperature where up to 120 states are 
kept in each renormalization step. 

First, we consider the Hamiltonian for 77 = (referred 
to as "isotropic" point below) where Eq. ([!]) simplifies to 
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FIG. 1: The main figure (inset) shows the temperature depen- 
dence of nearest neighbor correlations (Si ■ Si+i) ((t* ■ Tj+i)) 
for 77 values 0.0, 0.07, 0.08, 0.1, 0.11, 0.12, 0.15. For com- 
parison the spin only case is shown where Jij is fixed by its 
zero temperature expectation value for r\ = 0.15. The circle 
(square) at T = in the inset denotes — In 2 + 1/4 (—3/8). 



ff = ^(S,-S, +1 + l) 
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Note, that the model has a SU(2)xSU(2) symmetry but 
not the higher SU(4) symmetry group as possible in the 
5 = 1/2 case 0j. Regarding only a single bond, the low- 
est energy — J/2 is obtained if the orbitals form a singlet 
with (t, • Ti + \) = — 3/4 and the spins are ferromagneti- 
cally aligned. The first excited state is given by a spin sin- 
glet/orbital triplet and is separated by a gap of J/4 only. 
Therefore a strong orbital dimerization in the ground 
state is expected leading to alternating ferromagnetic- 
antiferromagnetic spin exchange couplings, which would 
be given by —J/4 and +J/8 for perfect dimerization. 
Due to translational invariance the corresponding classi- 
cal ground state would be fourfold degenerate. We also 
might expect mixed spin-orbital excitations playing an 
important role in thermodynamics. In Fig. |l| the temper- 
ature dependence of the nearest neighbor spin-spin corre- 
lation function (Si ■ 5i+i) and orbital-orbital correlation 
function (n ■ Tj+i) for different r] values is shown. First, 
we want to discuss in more detail the case r\ — and turn 
to finite Hund's coupling later on. Note, that in the zero 
temperature limit the nearest neighbor orbital correla- 
tion approaches —3/8, the value for perfect dimerization. 
In the high temperature limit the orbitals are completely 
disordered and the spin pair correlation is therefore neg- 
ative. While lowering the temperature orbital singlets 
are formed and (Si ■ 5i+i) increases due to the ferromag- 
netic spin interaction within an orbital singlet. The curve 
reaches a maximum at T/J 0.13 corresponding to the 
energy scale of the orbital gap (see below), and antifer- 
romagnetic quantum fluctuations of spins in neighboring 
orbital singlets are strongly enhanced and dominate be- 
low. This emphasizes the differences between a quantum 



model with 5=1 and the classical limit 5^1 in- 
vestigated by Shen et al. (T^| where fluctuations in the 
antiferromagnetic bonds completely vanish. 

An alternating ferromagnetic-antiferromagnetic spin 
exchange should also show up in the spin-spin corre- 
lation function (SqS*) showing 4 site periodicity (i.e. 
7r/2-oscillations) and in a spin dimer correlation function 
(SQSfS*S* +1 ) with long range order and 7r-oscillations. 
In Fig. @ numerical results for some leading CLs and 
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FIG. 2: The main figure shows the leading CLs times temper- 
ature in the isotropic case; the inset shows the corresponding 
wavevectors (in units of n) with circles denoting k — ±n/2. 

the corresponding wavevectors are shown. Because it 
is known from conformal field theory that CLs £ be- 
longing to critical excitations diverge as £ ~ 1/T in the 
low-tempcrature limit, CL times temperature versus tem- 
perature is plotted, thus separating correlation functions 
with long range order at T = from short range corre- 
lations. A dimer CL showing k — 7r oscillations diverges 
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as £ ~ exp(A D /T)/VT with A D w 0.08 in the low- 
temperature limit indicating long-range dimer order at 
zero temperature. Second largest at low temperatures is 
a spin-spin CL showing incommensurate oscillations ap- 
proaching 7r/2 at zero temperature. Thus, the numerical 
results are consistent with the picture above. Note also, 
that at high temperatures the orbitals have a short range 
antifcrromagnetic order and that a spin-orbital CL shows 
up, indicating the importance of mixed excitations. 

For a dimerized orbital pseudospin- 1 /2 as well as for 
an S — 1 ferromagnetic-antiferromagnetic alternating 
Heisenberg chain a gap in the excitation spectrum is 
expected as visible in the susceptibility data shown in 
Fig. [|a. The quadratic dispersion of a gapped sys- 
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FIG. 3: (a) Spin and orbital susceptibility for r\ = with low- 
temperature fits as specified in the text, (b) Specific heat for 
the isotropic and the finite Hund's coupling case . 

tern leads in ID to the low-temperature asymptotics 
X ~ exp(—A/T)/y/T. Using this function for a fit, 
a spin gap A s = (0.041 ± 0.002) J and an orbital gap 
A orb = (0.120±0.008)Jis obtained. Note that the orbital 
gap is much smaller than the singlet-triplet splitting J of 
an isolated orbital dimer. The reason is that in the cou- 
pled spin-orbital system the spins form a singlet if the or- 
bitals are excited to a triplet configuration and thus, the 
much smaller gap of the orbital triplet/spin singlet ex- 
citation appears in the orbital susceptibility data. More 



precisely, the thermal gap visible here is half of the spec- 
tral gap |^(|, so that we obtain A^ ctr ' ~ 0.24 J. This is 
very close to the value J/ 4 for an isolated bond, showing 
that the orbitals are nearly perfect dimers in the ground 
state. Dynamical coupling of two sectors affects also the 
spin gap value. In a simple picture of spin-2 objects cou- 
pled antiferromagnetically with J/32 jl4| , one would ex- 
pect the S = 2 Haldane gap of about 0.04(J/32) = J/800 
p7[ . The observed spin gap is in fact much larger due 
to the coupling between spins and orbitals: A spin exci- 
tation introduces additional antiferromagnetic couplings 
between the orbital singlets, thus raising also the energy 
of the orbital sector without destroying the dimer state. 

Now we want to discuss the effects of finite Hund's 
coupling. As easily understood from Eqs. (]]],§) an addi- 
tional, orbital free term R tends to stabilize a ferromag- 
netic spin order at large values of r\. Numerical data for 
i] = 0.15 indeed show that (Si ■ Si+i) is monotonously in- 
creasing and reaches I , the value expected for an uniform 
ferrochain, in the zero temperature limit (Fig. |l|). This 
leads to an exact cancellation of the orbital anisotropy 
terms in the ground state (see Eqs. |l]j|) and therefore to 
an isotropic antiferromagnetic coupling of the orbitals. 
From the Bethe ansatz result for the Heisenberg chain 
one expects (Tj ■ t,-+i) = — In 2 + 1/4 consistent with the 
numerical result (see Fig. [l]). Low temperature asymp- 
totics of the spin and orbital expectation values suggest a 
first order phase transition between two possible ground 
states: Spin/orbital dimer phase and spin saturated fer- 
rophase. A critical coupling rj c ~ 0.11 separating them 
at zero temperature is just slightly below realistic values 
for vanadium oxides ~ 0.12 Q. 

The specific heat (Fig. ||b) shows a two peak struc- 
ture. The high temperature hump is due to orbital and 
mixed spin-orbital excitations whereas the low-T peak 
is of spin origin. With increasing Hund's coupling, the 
hump is shifted to higher temperatures, whereas the low- 
T peak is smaller and shifted to lower temperatures. It 
is striking that the high-T hump contains such a large 
entropy weight. If it is due to orbital excitations only 
its entropy weight would be In 2 but in fact it is approx- 
imately 1 . This could be explained by a formation of an 
intermediate dimer state at finite temperature. In this 
case the entropy weight of the hump is given by the to- 
tal entropy In 6 minus a spin entropy of S — 2 dimers 
given by 1/2 In 5 consistent with the numerical data. To 
support this scenario of a temperature driven dimeriza- 
tion further, we have calculated again the leading CLs 
(Fig. |i|). The dimer CL with 7T-oscillations is indeed the 
largest CL in a certain temperature range and hence, 
dimerization is in this temperature range the leading in- 
stability towards the onset of long range order in a 3D 
system. Note also, a large orbital CL with 7T-oscillations 
indicating antiferromagnetic orbital order. At low tem- 
peratures there is a non-oscillating spin CL present as 
expected for an uniform ferrochain, however, this ferro- 
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FIG. 4: The leading CLs times temperature for r\ — 0.15. 
Note, that all shown CLs have commensurate oscillations over 
the entire temperature range. 



magnetic CL is strongly suppressed at finite temperature 
due to a disorder effect of low- lying dimerized states. The 
competition between uniform and dimerized states is ev- 
ident also from Fig. In comparison to the spin only 
case, where Jjj is fixed by its zero temperature expec- 
tation value for r/ = 0.15 and therefore all dimerization 
effects are omitted, (Si ■ S'i+i) shows an unusual rapid 
decay of local spin correlations with temperature. 

Finally, we address the relevance of these findings to 
the intermediate phase of YVO3. Obviously we have to 
take the coupling between different c-axis chains into ac- 
count. The quasi ID character of orbital physics is sup- 
ported by the C-type spin structure leading to a weak 
interchain orbital coupling proportional to 77 only [ fl2| . 
Interchain spin couplings can be explicitly included in 
the numerical calculations as a mean field, AJ a i,(S z ). An 
estimation for J a b can be obtained in the following way: 
Within linear response the susceptibility of the 3D sys- 
tem is given by \3D = Xid/(1 - ^JabXm)- Taking 
Tni/J ~ 0.25 from experiment (a representative value 
of J ~ 40meV is used), demanding that X3D diverges at 
this temperature and using the numerical data for \id 
calculated at 77 = 0.15, we find J a b/J ~ 0.075 in good 
agreement with the experimentally observed value [ fl9|| . 
In Fig. [| the dimer correlation length with this inter- 
chain molecular field included is displayed, showing that 
the results above are only weakly affected. The temper- 
ature range where dimerization is the leading instability 
approximately coincides with the C-phase of YVO3. We 
therefore believe that the model explains why in YVO3 a 
C-type magnetic order is stable as an intermediate phase. 
The gain of entropy due to the dimerization lowers the 
free energy F = (H) — TS at finite temperature and 
stabilizes a formation of alternating weak and strong fer- 
robonds along one cubic axis (e.g. caxis). In other words, 
a dimerization of the ferromagnetic spin chain occurs due 
to the orbital Peierls effect, in which thermal spin fluc- 



tuations play a role of lattice degrees of freedom. The 
effect is dynamical in the present ID case. However, as 
discussed before a weak but finite coupling between dif- 
ferent c-chains is present in YVO3 which does not destroy 
the quasi ID spin-orbital physics but leads to the onset 
of long range dimer order. Doubling of the magnetic 
unit cell along the c axis is then expected. In fact, a 
modulation of ferrocouplings and related splitting of the 
magnon branches has recently been observed by Ulrich 
et al. |l9| in the C-phase of YVO3. As the dimerization 
is of electronic origin, concomitant lattice distortions are 
expected to be small because t^g orbitals are not bond 
oriented and only weakly coupled to the lattice. 

Summarizing, the one-dimensional spin-orbital model 
with S = 1 shows an intrinsic tendency towards dimer- 
ization. Large Hund's coupling can stabilize a uniform 
ferromagnetic ground state but strong dimer correlations 
develop again at finite temperature due to the orbital 
Peierls effect. Dimer states modulate spin and orbital 
bonds and release high entropy. The quasi ID physics 
is not destroyed when realistic interchain couplings are 
taken into account on a mean field level. The temper- 
ature driven dimerization in the ID spin-orbital model 
could explain the stability of the intermediate C-type 
phase in YVO3. 
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lows: The elementary excitations are orbital "spinons" 
each with a gap A, but they always appear in pairs so 
that A^ c r ' = 2A. In the thermodynamic data the ac- 
tivated behaviour is characterized by the energy A of an 



elementary excitation irrespective of whether these exci- 
tations appear in pairs or not. 



